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LONG-TERM  GOALS 

The  overall  technical  goal  is  to  implement  a  1/25°  horizontal  resolution  global  ocean  prediction  system 
based  on  the  HYbrid  Coordinate  Ocean  Model  (HYCOM)  with  tides  and  dynamic  sea  ice.  The 
scientific  goals  include  but  are  not  limited  to  a)  evaluation  of  the  internal  tides  representation  in 
support  of  field  programs,  b)  data  assimilation  in  the  presence  of  tides,  c)  evaluation  of  the  model’s 
ability  to  provide  useful  boundary  conditions  to  high  resolution  coastal  models,  d)  interaction  of  the 
open  ocean  with  ice,  e)  shelf-deep  ocean  interactions,  f)  upper  ocean  physics  including  mixed 
layer/sonic  depth  representation,  g)  mixing,  etc. 

OBJECTIVES 

Perform  the  R&D  necessary  to  develop,  evaluate,  and  investigate  the  dynamics  of  1/25°  global 
HYCOM  with  tides  coupled  to  CICE  (the  Los  Alamos  sea  ice  model)  with  atmospheric  forcing  only, 
with  data  assimilation  via  NCODA  (NRL  Coupled  Ocean  Data  Assimilation),  and  in  forecast  mode. 
Work  closely  with  NRL  Stennis  to  incorporate  advances  in  dynamics  and  physics  from  the  science 
community  into  the  HYCOM  code  maintained  by  the  Navy. 

APPROACH 

This  past  year,  a  series  of  HYCOM  configurations  is  used  to  (b)  investigate  the  interaction  of  the  open 
ocean  with  ice,  (b)  evaluate  the  AMOC  representation  and  (c)  evaluate  its  performance  under  the 
CORE  forcing.  The  current  HYCOM  development  is  the  result  of  collaborative  efforts  among  the 
Florida  State  University,  University  of  Michigan,  and  the  Naval  Research  Laboratory  (NRL)  as  part  of 
the  multi-institutional  HYCOM  Consortium  for  Data- Assimilative  Ocean  Modeling  (Bleck,  2002; 
Chassignet  et  al.,  2003;  Halliwell,  2004).  This  effort  was  initially  funded  by  the  National  Ocean 
Partnership  Program  (NOPP)  to  develop  and  evaluate  a  data-assimilative  hybrid  isopycnal-sigma- 
pressure  (generalized)  coordinate  ocean  model  (Chassignet  et  al.,  2009).  HYCOM  has  been  configured 
globally  and  on  basin  scales  at  up  to  1/25°  (~3.5  km  mid-latitude)  resolution.  More  details  on  the  latest 
simulations  can  be  found  at  http://www.hycom.org  and  in  the  separate  ONR  reports  on  NRL  activities 
by  A.  Wallcraft  and  the  U.  of  Michigan  activities  by  B.  Arbic. 
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RESULTS 


Arctic  Ocean 


The  Arctic  Ocean  model  runs  were  performed  with  the  HYbrid  Coordinate  Ocean  Model  (HYCOM) 
coupled  with  Los  Alamos  National  Laboratory  Community  Ice  Code  (CICE).  Two  model 
configurations  of  the  Arctic  Ocean  were  used  in  the  study:  model  with  0.08°  (ARCcO.08)  and  with 
0.72°  (ARCcO.72)  horizontal  spacing  (Fig.l).  The  vertical  grid  has  32  hybrid  layers. 
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Fig.  1.  Model  domain  and  horizontal  resolution  (km)  of  ARCcO.08. 

High-latitude  storms  and  water  mass  formation  processes  in  the  Nordic  Seas 


The  Nordic  Seas  (Greenland,  Norwegian,  Iceland,  and  Barents  Seas)  is  a  key  region  for  maintenance 
of  the  Arctic  Ocean  thermohaline  structure.  This  is  the  region  of  intense  water  mass  formation  through 
cooling,  brine  rejection,  and  mixing  of  Arctic  Ocean  and  Atlantic  waters.  Misrepresentation  of  these 
processes  in  the  ocean  models  will  eventually  result  in  biases  of  thermohaline  structure  of  the  Arctic 
Ocean.  Ocean  processes  involved  in  the  water  mass  formation  in  that  area  are  largely  controlled  by 
surface  winds.  The  most  energetic  winds  in  the  Nordic  Seas  are  associated  with  storms  or  cyclones. 
There  are  two  major  categories  of  cyclones  in  that  region:  the  large-scale  and  meso-scale  low-pressure 
systems.  Meso-  and  small-scale  cyclones  are  low-pressure  systems  with  spatial  scales  in  O(<10  km) 
and  time  scales  from  several  hours  to  days.  A  subtype  of  this  class  of  cyclones  is  Polar  Low.  These 
lows  are  very  intense  maritime  low-pressure  system  with  strong  near-surface  winds.  Polar  Lows  form 
over  the  sea  and  predominantly  during  the  winter  months.  The  small-scale  cyclones  are  difficult  to 
observe  due  to  their  short  life  cycle  and  small  size.  It  is  speculated  that  the  role  of  the  meso-scale  low- 
pressure  systems  in  ocean  processes  is  as  large  as  that  of  the  large-scale  cyclones,  mostly  due  to  their 
higher  frequency  and  stronger  winds.  Meso-  and  small-scale  cyclones  over  the  Arctic  Ocean  are  poorly 
represented  in  the  available  observational  reanalysis  data.  The  following  wind  products  are  analyzed  in 
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this  study:  Cross-Calibrated  Multi-Platform  surface  wind  data  (CCMP);  the  NCEP/NCAR  Reanalysis  2 
(NCEPR);  NCEP  Climate  Forecast  System  Reanalysis  (CFSR);  and  Arctic  System  Reanalysis  (ASR). 
Analysis  of  the  wind  products  revealed  that  meso-scale  features  were  not  represented  in  either  of  the 
wind.  Even  large-scale  storms  have  different  representation  in  the  wind  products  (Fig.  3).  A  suite  of 
model  experiments  forced  by  different  winds  was  conducted  for  2005-2006. 


Fig.  3.  Different  representation  of  the  large-scale  cyclone  over  the  Greenland  Sea  in  the  wind 
products.  Note  the  difference  in  the  maximum  winds. 


We  find  that 

1 .  Different  wind  products  have  obviously  different  representation  of  storms  in  high  latitudes 
(Fig.  3).  This  impacts  all  wind-driven  process  in  the  Arctic. 

2.  The  upper-ocean  currents  from  the  model  experiments  forced  by  different  winds  differ  in 
magnitude  and  direction  (Fig.  4). 

3.  The  air-sea  fluxes  are  substantially  different  under  the  storms  in  the  simulations  driven  by 
different  winds  (Fig.  5). 

4.  Discrepancies  in  the  momentum  and  heat  fluxes  provided  by  different  winds  to  the  ocean 
impact  the  water  mass  formation  processes  in  the  model  (Fig.  6). 


3 


Fig.  4.  Winter  mean  upper-ocean  circulation  in  the  Nordic  Seas  from  the  ARCcO.08  simulation. 


Fig.  5.  Instantaneous  net  surface  heat  fluxes  (W/m2)  during  storm  shown  in  Fig.  3  in  the 
simulations  forced  by  winds  from  different  products  (indicated  in  the  titles). 


CFSR  NCEPR  CCMP  ASR 

Experiments 


Fig.  6.  Integrated  dense  water  production  (black  bars)  and  export  (gray)  out  of  the  Barents  Sea 
(left)  and  Greenland  Sea  (right)  in  the  simulations  driven  by  different  winds  (horizontal  axis). 


4 


Freshwater  pathways  in  the  Nordic  Seas 

The  Nordic  Seas  are  a  region  of  deep  water  formation.  Deep  waters  formed  in  the  convection  regions 
constitute  a  component  of  the  Meridional  Overturning  Circulation.  Several  studies  relate  intensity  of 
deep  convection  in  the  Nordic  Seas  to  freshwater  content  of  the  East  Greenland  Current  arguing  that  a 
substantial  fraction  of  this  freshwater  is  supplied  to  convection  regions.  However,  the  actual  pathways 
of  freshwater  from  the  East  Greenland  Current  to  the  convection  regions  are  not  well  known.  Of 
particular  interest  is  the  fate  of  freshwater  flux  from  Greenland  which  has  substantially  increased  over 
the  last  decade.  Can  this  freshwater  propagate  to  deep  convection  regions?  This  study  employs 
ARCcO.08  and  ARCcO.72  to  study  freshwater  pathways  with  origins  in  the  East  Greenland  Current. 
Model  experiments  are  used  to  analyze  the  role  of  wind  forcing  in  the  upper  ocean  circulation  and 
propagation  of  freshwater  in  the  region. 

To  date,  a  control  run  (with  no  Greenland  rivers)  has  been  performed  for  2006  with  (30-day)  and 
“without”  (4-year  relaxation  time  scale)  sea  surface  salinity  relaxation.  A  second  set  of  experiments 
has  be  run  with  monthly  high-resolution  Greenland  runoff  incorporated  into  the  model.  The  runoff 
data  are  provided  by  J.  Bamber  for  1999  -  present  time  period.  Test  experiments  with  and  without  SSS 
relaxation  have  been  conducted  for  2006.  Model  results  are  under  analysis.  It  is  noteworthy  that  for  the 
experiments  with  Greenland  Runoff,  changes  have  been  done  to  the  HYCOM  code  to  allow  river  water 
to  be  distributed  under  the  sea  ice. 


ARCcO.08-01 1  2006/08/1 5;  S  Z=-2m  ARCcO.08-051  2006/08/1 5;  S  Z=-2m 


Fig.  7.  Snapshots  of  the  upper  ocean  salinity  field  around  Greenland  on  August  15,  2006 from  the 
ARCcO.08  simulation  without  (left)  and  with  (right)  Greenland  Runoff  The  simulation  is  with  4- 

year  relaxation  time  scale. 

Relation  of  current  Arctic  climate  change  to  freshwater  flux  to  the  Nordic  Seas 

One  of  the  most  intriguing  manifestations  of  observed  climate  change  is  the  cessation  of  decadal 
Arctic  climate  variability  in  the  21st  century.  Between  1950  and  1996,  the  Arctic  atmospheric  and 
oceanic  circulation  alternated  between  anticyclonic  and  cyclonic  circulation  regimes  at  5  to  7  year 
intervals.  Since  1997,  however,  the  Arctic  system  has  been  dominated  by  a  16-year  anticyclonic 
regime  with  environmental  parameters  that  are  atypical  for  anticyclonic  climate  state.  In  the  20th 
century,  freshwater  and  heat  exchange  between  the  Arctic  Ocean  and  the  Nordic  Seas  were  self- 
regulated  and  their  interactions  were  realized  via  decadal  climate  auto-oscillations.  It  is  hypothesized 
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that  in  the  21st  century,  these  regular  oscillations  have  been  interrupted  as  a  result  of  additional 
freshwater  fluxes  from  Greenland.  The  excess  freshwater  from  Greenland  reduces  deep  convection  in 
the  northern  seas,  resulting  in  the  cessation  of  decadal  Arctic  climate  oscillations.  Analyzing  NAO 
index  and  air  temperature  (65-90°N)  time  series  we  have  found  that  decadal  changes  in  the  Arctic 
were  absent  during  Arctic  warming  in  the  1920s  -  1950s  when  Greenland  ice  melting  was  increased. 
Model  experiments  with  different  freshwater  fluxes  to  the  Nordic  Seas  (described  above)  are  currently 
used  to  investigate  this  hypothesis. 

AMOC  Variability 

The  state-of-the-art  approach  to  measuring  the  Atlantic  Meridional  Overturning  Circulation  (AMOC) 
is  to  continuously  observe  the  currents  and  hydrographic  fields  with  a  long-term,  moored  instruments 
array,  to  obtain  the  temporal  mean  value  of  the  AMOC  as  well  as  the  information  on  its  time  scale  and 
range  of  the  variability.  For  this  purpose,  the  joint  U.K.-U.S.  Rapid  Climate  Change-MOC  program 
(RAPID,  hereafter)  have  instrumented  a  complete  transbasin  section  along  26.5“N.  Daily  mean  AMOC 
transport  and  its  vertical  structure  have  been  derived  from  the  RAPID  data  since  April  2004  (e.g., 
Cunningham  et  al.,  2007;  Kanzow  et  al.,  2007,  2010;  Johns  et  al.,  2008,  2011;  Rayner  et  al.,  2011; 
McCarthy  et  al.,  2012).  The  RAPID  results  show  that  the  AMOC  is  not  slowly  varying,  as  previously 
thought,  but  rather  has  a  high  variability  on  intraseasonal  to  seasonal  time  scales.  In  addition,  there  is 
significant  variability  on  interannual  and  longer  time  scales.  One  question  that  arises  is  as  a  to  whether, 
or  to  what  extent,  the  variability  is  coherent  across  different  latitudes,  so  that  the  variability  observed 
at  one  latitude  is  representative  of  a  larger  scale  change  (Willis,  2010).  This  question  can  be  addressed 
using  a  high  horizontal  ocean  model  as  long  the  model  adequately  represents  the  basic  features  of 
AMOC  variability. 

Results  from  two  eddy-resolving  HYCOM  numerical  simulations  are  presented  here.  The  two 
simulations,  one  Atlantic,  the  other  global,  have  several  key  features  in  common.  Both  have  a 
horizontal  resolution  of  1/12'’.  Both  were  initialized  using  potential  temperature  (T)  and  salinity  (S) 
from  the  ocean  climatology,  the  Generalized  Digital  Environmental  Model  [GDEM,  Carnes,  2009], 
and  were  spun-up  from  rest  using  the  climatological  surface  forcing  from  the  ECMWF  reanalysis 
ERA40  [Uppala  et  al.,  2005].  After  spin-up  (10  years  for  the  global  model  and  15  years  for  the 
Atlantic),  both  use  the  3-hourly,  0.5°  Navy  Operational  Global  Atmospheric  Prediction  System 
(NOGAPS,  Rosmond  et  al.,  2002).  The  global  simulation  covers  2003-2012,  the  Atlantic  simulation 
covers  2004-2012.  As  a  subset  of  the  global  model,  the  Atlantic  model  domain  extends  meridionally 
from  28°S  to  the  Fram  Strait  at  80°N.  No  inow  or  outow  is  prescribed  at  the  northern  and  southern 
boundaries.  Within  a  buffer  zone  of  about  3°  from  the  northern  and  southern  boundaries,  the  3-D 
model  T,  S,  and  depth  of  isopycnal  interface  are  restored  to  the  monthly  GDEM  with  an  e-folding  time 
of  5-60  days  that  increases  with  distance  from  the  boundary.  Vertically,  the  Atlantic  simulation  has  the 
conventional  32  layers  whereas  the  global  model  has  41  layers.  The  additional  nine  layers  are  all  near 
the  surface  to  resolve  the  seasonal  thermocline  in  other  parts  of  the  global  ocean.  The  Atlantic  model, 
with  dfferent  reference  densities  near  bottom,  uses  some  of  the  nonexistent  bottom  layers  in  the  global 
model  and  has  a  finer  vertical  resolution  in  the  deep  Atlantic  Ocean.  This  change,  as  documented  in 
Xu  et  al.  (2012),  improves  the  vertical  structure  of  the  AMOC.  The  global  simulation,  which  includes 
a  coupled  sea-ice  model,  instead  of  the  simple  energy  loan  sea-ice  parameterization  in  the  Atlantic 
simulation,  shows  improved  water  property  (salinity  in  particular)  in  the  upper  layers  of  the  Labrador 
Sea.  However,  the  southward  flow  of  NADW  in  the  global  simulation  has  significantly  less 
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contribution  from  the  Nordic  Seas  overow  water  (NSOW),  probably  because  of  inadequate 
representation  of  the  model  water  property  in  the  Nordic  Seas. 

We  employed  an  adaptive  and  temporally  local  decomposition  method  to  separate  the  transport  time 
series  into  variability  modes  of  different  time  scale.  The  results  show  that  at  26.5°N,  our  simulations 
yield  most  of  the  observed  AMOC  variability  on  intraseasonal,  seasonal,  and  interannual  time  scales. 
The  variability  of  the  Florida  Current  transport  in  the  Atlantic  model  also  agrees  resonably  well  with 
the  observations.  Near  41‘’N,  our  model  results  represent  well  the  observed  variability  on  seasonal  time 
scales,  but  is  less  successful  on  interannual  time  scales.  On  a  long  time  scale,  the  observed  AMOC 
transports  at  26.5°N  declined  about  3  Sv  from  2004  to  2012.  The  Atlantic  model  does  not  simulate  this 
long-term  change,  but  the  global  model  does.  Thus,  the  change  may  represent  large-scale  dynamics 
beyond  the  Atlantic  model  domain.  Overall,  the  variability  of  the  monthly  mean  AMOC  transport  in 
both  the  global  and  the  Atlantic  simulations  is  lower  than  the  observations.  The  variability  of  the 
AMOC  transport  in  the  two  simulations  is  very  similar,  despite  differences  in  the  model  domain  and 
the  vertical  structure  of  the  southward  limb.  This  indicates  that  most  of  the  interannual  and  shorter 
term  variability  in  our  models  is  not  induced  by  dense  water  formation  to  the  north,  or  by  the  dynamics 
to  the  south  of  Atlantic  model  domain.  The  variability  more  likely  is  induced  via  the  same  atmospheric 
forcing  used  in  two  simulations,  particularly  the  wind.  This  is  in  accord  with  recent  works  of  Zhao  and 
Johns  (manuscript  submitted  to  JPO;  2013  AMOC  meeting  presentation)  and  Yang  (2013  AMOC 
meeting  poster).  They  show  that,  in  idealized  two-layer  models  that  contain  no  time  mean  AMOC, 
seasonal  and  interannual  variability  of  the  AMOC  transport  is  induced  by  variation  of  the  wind-driven 
subtropical  gyre.  Also,  Ekman  transports  directly  contribute  to  the  AMOC  variability  on  various  time 
scales,  such  as  the  abnormally  low  AMOC  transport  in  Winter  2009/2010  caused  by  the  very  low 
North  Atlantic  Oscillation  (Xu  et  ak,  2012;  Hakkinen  and  Rhines,  2013).  The  notion  of  wind  as  the 
primary  source  of  AMOC  varability,  at  least  within  interannual  time  scales,  highlights  the  need  for 
better  quantifying  of  the  wind  field  and  better  understanding  of  the  wind-driven  ocean  responses. 

When  AMOC  variability  at  different  latitudes  is  considered,  the  results  show  little  meridional 
coherence  on  an  intraseasonal  time  scale,  although  large  intraseasonal  variations  sometimes  extend 
over  long  distances  meridionally.  The  seasonal  variability  of  the  AMOC  is  connected  from  the 
subtropical  to  the  subpolar  North  Atlantic,  but  is  clearly  separated  from  the  subtropics  to  the  tropics. 

For  the  well-documented  low  transport  at  26.5°N  during  the  winter  of  2009/2010,  there  is  an  opposite 
response  in  the  subpolar  North  Atlantic.  All  the  northward  AMOC  transports  in  the  North  Atlantic 
come  from  the  South  Atlantic.  The  model  results  show  that  the  AMOC  varies  dramatically  on  seasonal 
time  scales  when  crossing  the  tropical  Atlantic.  This  seasonal  variation  is  primarily  due  to  Ekman 
transport  near  the  surface,  but  is  compensated  to  a  small  part  by  contribution  of  AAIW  below  the 
thermocline  (Fig.  8).  The  importance  of  the  Ekman  transport  in  the  tropical  Atlantic  is  highlighted 
further  in  heat  transport:  At  10°N,  the  Ekman  transport  explains  nearly  all  the  northward  heat  flux. 
Horizontally,  the  circulation  pattern  in  the  upper  layer  of  the  tropics  also  change  significantly  in 
association  with  the  seasonal  AMOC  variation.  Particularly  interesting,  in  August  the  NECC  gains 
transport  from  the  north  as  it  ows  eastward  and  recirculates  south  and  westward  to  form  a  semiclosed 
gyre.  This  pattern  differs  from  the  classical  picture  of  near-surface  circulation  in  the  tropical  Atlantic 
by  Richardson  and  Walsh  (1986),  but  explains  the  slightly  southward  net  transport  in  August.  These 
results  emphasize  the  value  of  eddy-resolving  simulations  in  providind  a  detailed  spatial  structure  of 
the  circulation  in  addition  to  transbasin  volume  and  heat  transport  estimates. 
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Fig.  8:  Monthly  mean  transports  of  volume  (a,  c)  and  temperature  (b,  d)  across  lO^N  as  a  function 
of  model  layers.  Circles  in  panel  (c)  denote  the  Ekman  transport  determined  from  model  wind 
stress.  Results  are  based  on  a  9-year  mean  of  a  1/1 2° Atlantic  simulation  in  2004-2012. 

CORE  experiments 

The  aim  of  the  Coordinated  Ocean-ice  Reference  Experiments  (CORE)  is  to  compare  results  from 
several  ocean-ice  model  using  the  same  atmospheric  forcing  to  create  a  basis  that  could  be  used  to 
provide  a  mechanistic  understanding  of  the  role  of  the  ocean  in  the  climate  system  and  to  better 
quantify  the  importance  of  numerical  choices.  HYCOM,  coupled  with  the  Los  Alamos  CICE,  was  set¬ 
up  following  the  CORE_v2  protocol.  The  protocol  consists  of: 

•  an  interannually  varying  atmospheric  forcing  (corresponding  to  the  period  1948-2007)  from 
Large  and  Yeager  (2009)  for  5  consecutive  cycles. 

•  an  SSS-relaxation  corresponding  to  4years/50m 

•  a  surface  salinity  global  normalization 


The  HYCOM-CICE  global  configuration  uses  version  2.2.74  of  HYCOM  and  version  4.0  of  CICE.  It 
has  a  500x382  points  tripolar  grid  of  ~  0.72°  resolution  and  32  sigma-2  hybrid  vertical  layers.  The  sea- 
ice  model  provides  the  ocean  model  with  ice-ocean  stress,  SST  and  ice  cover.  The  ocean-ice  fluxes  are 
calculated  using  the  energy  loan  ice  model  built  in  HYCOM.  The  initial  conditions  and  the  surface 
salinity  climatology  is  the  General  Digital  Environmental  Model  (GDEM4;  Carnes  et  ah,  2010)  and 
the  turbulent  air-sea  fluxes  are  computed  using  the  Kara  bulk  formulae  (Kara  et  ah,  2004).  To  account 
for  the  temperature  bias  generated  by  this  particular  bulk  formulation,  a  constant  flux  offset  of — lOW 
is  applied  everywhere.  To  evaluate  the  sensitivity  to  the  surface  salinity  relaxation,  two  experiments 
were  performed:  1)  we  apply  a  30days/50m  SSS  relaxation  (SSS30d)  and  2)  we  apply  4years/50m  SSS 
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relaxation  (SSS4y).  In  both  cases,  a  surface  salinity  global  normalization  is  also  applied  (i.e.  no  net  salt 
flux).  The  model  is  first  run  30  years  (1948-1978)  to  evaluate  each  experiment’s  performance. 


Temperature 


Fig.  9.  Evolution  of  the  global  temperature  (top)  and  maximum  of  the  AMOC  (Sv;  bottom)  at 
26.5‘>N for  SSSSOd  (red),  SSS4y  (blue),  CICE-SSSSOd  (green),  CICE-SSS4y  (dashed-red)  and  CICE- 

varyingSSS  (yellow). 

Both  experiments  exhibit  a  similar  increase  of  the  global  temperature  (~0.15°C,  Fig.  9  top)  and  SSH 
(~9cm,  not  shown)  after  30  years.  A  strong  fresh  bias  is  observed  in  SSS4y  in  the  southern  Ocean  and 
the  Arctic  Ocean  (Fig.  10,  right).  This  fresh  bias  in  the  Arctic  then  partly  feeds  the  subpolar  gyre  in 
SSS4y  resulting  in  a  continually  decreasing  Atlantic  Meridional  Overturning  Circulation  (AMOC) 
from  17  Sv  at  the  beginning  to  5  Sv  after  30  years  (Fig  9,  bottom).  In  SSS30d,  the  surface  salinity 
relaxation  prevents  these  polar  fresh  biases,  keeping  the  SSS  of  the  supolar  close  to  the  climatology 
(Fig.  10,  left).  Consequently,  while  we  still  do  observe  a  decrease  of  the  AMOC  in  SSS30d,  it 
stabilizes  close  to  ~9  Sv  after  10-15  years. 


I 
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Fig  10.  Surface  salinity  biases  after  30  years  for  SSS30d  (left)  and  SSS4y  (right). 
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A  large  part  of  the  bias  in  SSS4y  was  traced  to  the  re-calculation  of  the  ocean-ice  fluxes  based  on 
solely  the  SST  and  ice-cover  using  the  energy  loan  model  built  in  HYCOM.  We  therefore  performed 
another  set  of  experiments  to  evaluate  the  impact  of: 

•  the  ocean-ice  surface  fluxes  as  provided  by  CICE  instead  of  recalculating  them  in  HYCOM  via 
the  energy  loan  model 

•  the  Levitus  climatology  for  the  initial  conditions  and  SSS  relaxation  instead  of  the  GDEM4 
climatology  to  better  balance  the  initial  ocean  temperature  with  the  COREv2  atmospheric  heat 
fluxes  and  consequently  improve  the  global  temperature  drift. 


Total  ice  extent  {1e12  m2)  -South  Pole- 

25  ^  ^  ^ 


-  Flux  exchange  sss4y 

_ Observation  Mean  1 979-201 0 

12  2 

Fig  11:  Evolution  of  the  southern  pole  ice  extent  (10  m  )  from  NSIDC-SSMI  (obs.)  in  dashed 
black,  CICE-SSSSOd  in  red  and  CICE-SSS4y  in  black. 

First,  we  performed  a  new  60-year  simulation  keeping  the  SSS  relaxation  strong  at  30  days  (CICE- 
SSS30d),  but  with  the  ocean-ice  fluxes  exchanged  between  HYCOM  and  CICE  and  intialized  from  the 
Levitus  climatology.  Results  are  similar  to  SSS30d  with  a  slight  decrease  of  the  global  temperature 
trend  (+0.1 7°C  instead  of  +0.23°C;  Fig.  9).  Decreasing  the  SSS  relaxation  to  4years/50m  (CICE- 
SSS4y)  leads  to  a  salty  bias  over  the  poles  that  contributes  to  a  strong  AMOC  but  to  the  detriment  of 
ice  formation  at  the  southern  pole  after  15  years  (Fig  11).  A  spatially  varying  SSS  relaxation  was 
implemented  with  a  stronger  relaxation  (6  months/50m)  to  the  Southern  Ocean  region.  The  CICE- 
varyingSSS  experiment  is  able  to  sustain  an  ice  sheet  at  the  South  Pole  for  the  whole  60  years  of 
simulation  with  the  global  temperature  stabilizing  after  20  year  around  3.70°C  (Fig.  9  top).  The 
maximum  strength  of  the  AMOC  at  26.5°N  for  this  last  simulation  reaches  now  12.5  Sv  at  the  end  of 
the  E'  cycle  (Fig.  9  bottom),  within  the  range  of  most  CORE  simulations  (Danabasoglu  et  al.,  2013). 
CICE-varying  has  been  carried  out  for  the  full  length  of  the  CORE  protocol  (5  cycles)  and  shows  a 
stabilization  of  the  global  temperature  at  3.7  EC,  exhibits  a  positive  bias  of  salt  in  the  Arctic  (~lpsu),  a 
sustained  ice  sheet  on  both  Poles,  but  a  strong  variability  of  the  strength  of  the  AMOC.  However,  the 
CICE-varyingSSS  presents  a  maximum  AMOC  of  14  Sv  average  over  the  last  20  years  of  the  last 
cycle,  well  within  the  average  of  the  other  CORE -It  models  (Danabasoglu  et  al.,  2013).  The  model 
outputs  have  been  submitted  to  several  teams  of  the  CORE-II  project  and  we  are  working  on  several 
papers  to  be  published  in  peer-review  journals.  All  experiments  have  been  carried  out  on  the  NAVY 
computers  Einstein  and  Kilrain. 
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IMPACT/APPLICATIONS 


The  1/25°  (3.5  km  mid-latitude)  resolution,  first  used  in  some  FY09  global  HYCOM  simulations,  is 
the  highest  so  far  for  a  global  oeean  model  with  high  vertical  resolution.  A  global  ocean  prediction 
system,  based  on  1/25°  global  HYCOM  with  tides,  is  planned  for  real-time  operation  starting  in  2015. 
At  this  resolution,  a  global  ocean  prediction  system  can  directly  provide  boundary  conditions  to  nested 
relocatable  models  with  ~1  km  resolution  anywhere  in  the  world,  a  goal  for  operational  ocean 
prediction  at  NAVOCEANO.  Internal  tides  and  other  internal  waves  can  have  a  large  impact  on 
acoustic  propagation  and  transmission  loss  (Chin-Bing  et  ah,  2003,  Wam-Vamas  et  ah,  2003,  2007), 
which  in  turn  significantly  impacts  Navy  anti-submarine  warfare  and  surveillance  capabilities.  At 
present,  regional  and  coastal  models  often  include  tidal  forcing  but  internal  tides  are  not  included  in 
their  open  boundary  conditions.  By  including  tidal  forcing  and  assimilation  in  a  fully  3-D  global  ocean 
model  we  will  provide  an  internal  tide  capability  everywhere,  and  allow  nested  models  to  include 
internal  tides  at  their  open  boundaries. 
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